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Abstract 

We present a complete calculation of the electroweak 0(a s a B ) corrections to three-jet produc- 
tion and related event-shape observables at electron-positron colliders. The Z-boson resonance is 
described within the complex-mass scheme, rendering the calculation valid both in the resonance 
and off-shell regions. Higher-order initial-state radiation is included in the leading-logarithmic 
approximation. We properly account for the corrections to the total hadronic cross section and for 
the experimental photon isolation criteria. To this end we implement contributions of the quark- 
to-photon fragmentation function both in the slicing and subtraction formalism. The effects of the 
electroweak corrections on various event-shape distributions and on the three-jet rate are studied. 
They are typically at the few-per-cent level, and remnants of the radiative return are found even 
after inclusion of appropriate cuts. 
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1 Introduction 



Jet production in e + e _ annihilation provides an ideal environment for studies of the dynamics of the 
strong interaction, described by the theory of quantum chromodynamics (QCD) [1]. The kinematical 
distribution of jets closely reflects the parton-level kinematics of the event. Consequently, the first 
observation of three-jet final states at DESY PETRA [2], produced through quark-antiquark-gluon 
final states [3] , provided conclusive evidence for the existence of the gluon. Jets are defined through a 
jet algorithm, which is a procedure to recombine individual hadrons into jets using a distance measure, 
resolution criterion and recombination prescription. The theoretical description applies the same jet 
algorithm to partons in the final state. Closely related to jet cross sections are event-shape distributions. 
Event-shape variables measure certain geometrical properties of hadronic final states, and can equally 
be calculated in perturbative QCD from partonic final states. 

Jet cross sections and event-shape distributions were studied very extensively at e + c~ colliders [4], 
and high-precision data are available from the LEP experiments ALEPH [5], OPAL [6], L3 [7], DEL- 
PHI [8], from SLD [9] at SLAC and from JADE [10] at DESY PETRA. The theoretical description of 
these data in perturbative QCD contains only a single parameter: the strong coupling constant a s . By 
comparing experimental results with the theoretical description, one can thus perform a measurement 
of a s from jet cross sections and event shapes. For a long period, the perturbative description of these 
observables was based on next-to- leading order (NLO) [11-13] in perturbative QCD, improved by the 
resummation of next-to-leading-logarithmic corrections [14] to all orders. The uncertainty on these 
theoretical predictions from missing higher-order terms results in a theory error on the extraction of 
a s , which was quantified to be around five per cent, and thus larger than any source of experimental 
uncertainty. 

Owing to recent calculational progress, the QCD predictions for event shapes [15,16] and three-jet 
production [17, 18] are now accurate to next-to-next-to-leading order (NNLO, a 2 af) in QCD pertur- 
bation theory. Inclusion of these corrections results in an estimated residual uncertainty of the QCD 
prediction from missing higher orders at the level of well below five per cent for the event-shape dis- 
tributions, and around one per cent for the three-jet cross section. Using these results (combined [19] 
with the previously available resummed expressions), new determinations of a s from event-shape and 
jet production data were performed, resulting in a considerable improvement of the theory uncertainty 
to three per cent from event shapes [20] and below two per cent from jet rates [21]. A further im- 
provement can be anticipated for the event shapes from the resummation of subleading logarithmic 
corrections [22]. 

At this level of theoretical precision, higher-order electroweak effects could be of comparable mag- 
nitude. Until recently, only partial calculations of electroweak corrections to three-jet production and 
event shapes have been available [23], which can not be compared with experimental data directly. In 
a previous work [24], we briefly reported our results on the first calculation of the NLO electroweak 
(a 3 a s ) corrections to three-jet observables in e + e~ collisions including the quark-antiquark-photon 
(qq-y) final states. Here, we describe this calculation in detail and perform extensive numerical studies 
on the impact of the electroweak corrections to three-jet-like observables at different e + e~ collider 
energies. 

The full NLO electroweak (a 3 a s ) corrections to jet observables contain four types of contributions: 
genuine weak corrections from virtual exchanges, photonic corrections to quark-antiquark-gluon (qqg) 
final states, gluonic corrections to qq-f final states, and QCD/electroweak interference effects in qqqq 
final states of identical quark flavour. The latter were not included in our previous work, but turn out 
to be numerically negligible as anticipated. 

Any jet-like observable at NLO in the electroweak theory receives virtual one-loop corrections and 
contributions from real photon radiation. Experimental cuts on isolated hard photons in the final 
state allow to suppress these real photon contributions. However, photons radiated inside hadronic 
jets can often not be distinguished from hadrons (like neutral pions), and are thus not removed by 
experimental cuts. The real photon contribution at NLO thus results from a complicated interplay of 
jet reconstruction and photon isolation cuts. Through the isolated-photon veto, these observables are 
sensitive to final-state particle identification, and thus to fragmentation processes. In our case, we must 
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include a contribution from quark-to-photon fragmentation [25] to obtain a well-defined and infrared- 
safe observable. Since our calculation is among the very first to perform electroweak corrections to jet 
observables with realistic photon isolation cuts, we describe the relevant calculational aspects in detail 
below. 

To define the observables considered here, we describe in Section 2 the jet-clustering algorithms 
used in c + c~ annihilation and the standard set of event-shape variables. In this section, we also review 
the current description of these observables in perturbative QCD. The calculation of NLO electroweak 
corrections is outlined in Section 3, where we describe the calculation of the virtual and real corrections 
in detail. The real corrections contain infrared divergences from unresolved photon and gluon radiation. 
These infrared divergences cancel against similar divergences in the virtual corrections. To accomplish 
this cancellation, it is, however, necessary to extract them analytically from the real corrections, which 
is done by a slicing or subtraction procedure (described in Section 4) . Our results are implemented into 
a parton-level event generator (described in Section 5), which allows the simultaneous evaluation of 
all event-shape variables and jet cross sections. Numerical results for jet production and event-shape 
distributions for e + e~ collision energies at LEP1, LEP2 and a future linear collider are presented 
in Section 6. At energies above the Z peak, we observe non-trivial kinematical structures in the 
distributions. It is shown that these structures are a remnant of the radiative-return phenomenon, 
resulting from a complicated interplay of event-selection and photon isolation cuts applied in the 
experimental definition of the observables. Finally, we conclude with an outlook on the impact of 
electroweak effects on future precision QCD studies at c+c~ colliders in Section 7. 



2 Jet observables 

A commonly used method for reconstructing jets was originally introduced by the JADE group [26]. 
The algorithm is based on successive combinations. In a first step, each observed particle is listed as 
a jet. In the next step, a resolution parameter y^ is calculated for each particle pair, and the particle 
pair leading to the smallest value of yij is combined into a single pseudo-particle. This yields a new list 
of jets, and the algorithm proceeds with step two. The procedure is repeated until no pair of particles 
is left with < y C ut, where y cut is a preset cut-off. 

Different proposals exist in the literature in how to define y^ (see e.g. Ref. [27]). The original 
JADE definition reads 

2^^(l-cos%) 

VijJ = ^2 ' (. 2A > 

where Ei is the energy of the i-th particle, cos % the angle between the particles, and E v - ls the total 
visible hadronic energy in the event. Improving upon this definition, different jet resolution parameters 
have been proposed. Most widely used at LEP was the kr or Durham algorithm [28,29], which defines 

2min (^ 2 ,^)(l-cos%) 

y»j,D = ^2 • I 2 - 2 ) 

vis 

In addition to the choice of jet resolution parameter, there also exist different ways of combining the 
four- momenta of the two particles with the lowest y^ to one four- momentum p^. In the so-called 
^-scheme one simply adds the two four- momenta, leading to = Pi + Pj . In the P-scheme the 
invariant mass of the pseudo-particle is set to zero by rescaling the energy 

Pij =Pi+Pj, Eij = \pi+pj\. (2.3) 

In the P -scheme, (2.3) is used to construct the resulting four-momentum, however, after each recom- 
bination E v i s is recalculated. Finally, in the i? -scheme the three-momentum rather than the energy 
is rescaled. 

Since an event containing three jets is due to the emission of a gluon off an (anti-)quark at a large 
angle and with significant energy, the ratio of the number of observed three-jet to two-jet events is, in 
leading order, proportional to the strong coupling constant. In general, the n-jet rate R n (y), which 
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depends on the choice of the jet resolution parameter y = y cui , is defined through the respective cross 
sections for n > 2 jets 

R n (y,VS) = ^, (2.4) 



Chad 

such that 



n=l 

and yfs is the centre-of-mass (CM) energy. In order to characterise the topology of an event a large 
number of observables have been developed. Most of them require at least three momenta of final-state 
particles to be non-trivial. In the following we introduce six variables which have been extensively used 
in experimental analyses: thrust T [30,31], the normalised heavy-jet mass p [32], the wide and total jet 
broadenings B w and B T [33,34], the C-parameter [35,36], and the transition from three-jet to two-jet 
final-state using j/ ijiD [28,29,37,38]. 

• Thrust is defined through _ 

T = max ( 2 - 6 ) 

\Pi\ 

where pi is the three-momentum of the z-th particle, and n is varied to maximise the momentum 
flow in its direction, yielding the thrust axis. 

• Every event can be divided into two hemispheres Hi and H 2 by a plane perpendicular to the 
thrust axis. In each hemisphere Hi one can calculate the invariant mass Mf, the larger of which 
yields the heavy-jet mass 

M h 2 ad = max (Ml M 2 2 ) , (2.7) 

and the normalised heavy-jet mass 

M 2 

(2.8) 

vis 

• Using the definition of the hemispheres from above, one can calculate the hemisphere broadenings 

Y^itzM \Pj x n\ 

B <= 7zM ■ < 2 - 9 » 

The wide and total jet broadenings By/ and Bt are then obtained through 

B w = max(B 1 ,B 2 ), B T =B l +B 2 . (2.10) 

• Starting from the linearised momentum tensor 

e^^Eff. M-u.3, (in) 

and its three eigenvalues Ai, A 2 , A 3 , the C-parameter is defined through 

C = 3(A 1 A 2 + A 2 A 3 + A 3 A 1 ). (2.12) 

• The jet transition variable Y3 is defined as the value of the jet resolution parameter for which an 
event changes from a three-jet-like to a two-jet-like configuration. 

In the following we often denote event-shape observables generically y. Taking (1 — T) instead of T 
the two-particle configuration is located at y = for all the event-shape variables defined above. 
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Figure 1: Lowest-order diagrams for c+c — > qqg. 



2.1 Event shapes and jet rates in perturbation theory 

At leading order (LO), 0(a 2 a s ), the first process that occurs at tree level in e + e~ annihilation is 
gluon radiation off a quark or antiquark (see Fig. 1). As mentioned above, by comparing the measured 
three-jet rate and event-shape observables with theoretical predictions, one can determine a s . 

In perturbation theory up to next-to-next-to-leading order (NNLO) in QCD, the expansion of 
a distribution in the generic observable y at CM energy y/s for renormalisation scale /i = \fs an d 
a s = a s (s), normalised to the Born cross section ao(s) of the process e + c~ — > qq is given by 



Ida fa s \ dA ,/a s \ 2 dB ,/a s \ 3 dC . 4 v 

- dy--{^)dy- + {^) d^+fcJ -dy- + °^' (2 ' 13) 



where A, B 7 and C denote the QCD contributions of LO, next-to- leading order (NLO), and NNLO. 
The experimentally measured event-shape distribution is normalised to the total hadronic cross section 
(Thad, which for massless quarks reads 

-had = <ro (l + (g) K X + (g) 2 K 2 + O (a 3 )) , (2.14) 
1 da /a s \ dA /a s \ 2 dB /a s \ 3 dC' 



such that 



/a s \dA /a s \ z d£> fa s \- i dU , N .„ 

CThad dy y ^ + y ^ + y ^ + ° ™ - ( 2 - 15 ^ 



where 



A = A, B = B — AK\ , C = C-S^i + AA"? - A^. (2.16) 

The coefficients in (2.16) up to NLO have been calculated in Refs. [11-13]. Furthermore, kinematical- 
ly-dominant leading and next-to- leading logarithms have been resummed [14,39], and non-perturbative 
models of power-suppressed hadronisation effects have been included [40-43] to increase the theoretical 
accuracy. Recently the first NNLO calculations have been completed [15-18], and the matching of next- 
to-leading logarithms and next-to-next-to-leading logarithms to the fixed-order NNLO calculation has 
been performed [19,22]. These results have subsequently been used in precision determinations [20-22] 
of the strong coupling constant a s . 

With regard to jet rates, fixed-order calculations are known up to next-to-next-to-next-to-leading 
order (N 3 LO) in QCD for the two-jet rate [17,44-46], up to NNLO for the three-jet rate [11-13,17,18], 
and up to NLO for the four-jet rate [47-51]. 

NLO electroweak (EW) corrections could be of comparable magnitude as the NNLO QCD cor- 
rections and are therefore worth further consideration. The factorisable EW corrections have been 
calculated in Ref. [52] and a further step towards the full NLO EW corrections has been made in 
Ref. [23]. In this work we describe the first calculation of the complete NLO EW corrections to the 
normalised event-shape distributions. First results of this calculation on the thrust distribution and 
the three-jet rate at y/s — Mz have been presented in Ref. [24]. 

In analogy to the QCD corrections, we write the total hadronic cross section including 0(a) cor- 
rections as 

CThad = a o (l+(£)<5 CT ,i+0(a 2 )), (2.17) 
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and the expansion of the observable da/dy as 



(Tq dy V27T/ dy V27T/ dy \2ttJ \2tt J dy 



(2.18) 



where the LO purely electromagnetic contribution (5 7 arises from tree-level quark-antiquark-photon 
(qq-y) final states without a gluon* and 5a comprises the NLO EW corrections to the distribution 
da/dy. 

Normalising (2.18) to CThad yields 



1 da 



V2tt 



cL4 
dy 



CThad dy 

Hence, the full O(a) EW corrections' 1 ' are given by 
dfew 



+(3t+s)e)(t-^) +o( " 2 » 



dy 



/ a \ d<5 7 / a \ / a s \ ( d$A dA \ 
V2W ~dy~ + Vh^J Vhi) [~dy~ ~ ~dy~ aA ) ' 



(2.19) 



(2.20) 



In order to obtain a sensible ratio, all three contributions have to be evaluated using the same event- 
selection cuts. 

The EW corrections to both <7h a d an d the distribution in y contain large corrections due to initial- 
state radiation (ISR). Since these are universal, they partially cancel in the third term in (2.19), leaving 
only a small remainder. If we include higher-order leading-logarithmic (LL) ISR effects in both CThad 
and the distribution in y, this leads to 



Chad = a 



( 1 + (S)^ + (s) 2 ^ + °(« 2 ))' 



and 



1 da f a s\ dA / a \ d<5 7 / a \ / a s \ d$A ( a \ 2 / a s 
ctq dy _ \~2tt) du \2tt) ~~dv \~2~k) \2tt) ~~dv~ \~2tt) \2^ 



dA 
^2irJ dy 



dS 
dy 



dSA 
dy 



dS 



(2.21) 
(2.22) 

where <5 ct .>2,ll and <5a,>2,ll contain leading-logarithmic (LL) terms proportional to a n \n n (s/m 2 ) with 
n > 2, as defined in Section 4.4. Here O (a 2 ) stands for two-loop electroweak effects without the 
enhancement of leading ISR logarithms. For the normalised distribution this results in 



A,>2,LL 

dy 



1 da 
CThad dy 



/ as \ dA fa_\ d^ fa_\ / \ / d^ _ cL4 \ 
V2tJ dy V2tt/ dy V2tt/ \2ir) \ dy dy aA ) 



dA fa 
~d~y~^\2Tx 



dS A . 



>2.LL 



dy 



dA s 

-^<V,>2,LL 

dy " 



dS A . 



l.LL 



dy 



dy 



(2.23) 



where arid <j ct .i.ll denote the LL contributions contained in the NLO results. Their second- 

order effect results from taking the ratio of two NLO-corrected quantities. Therefore, the higher-order 
LL corrections read 



dS 



EW,LL 

dy 



-(s) (s) 



dS 



A,>2.LL d^4 
- 7-0a. 



dy 



dy 



>2,LL 



dA 

dy 



°cr,l,LL 



dS 



A,1,LL 

dy 



LL 



(2.24) 



Due to the universality of ISR, the terms in the first and in the second parenthesis in (2.24) separately 
cancel each other numerically to a large extent. 

The same decomposition as applied here for event-shape distributions holds also for the three-jet 
rate, normalised to CThad- 



* Since the event-shape observables are calculated from parton momenta, the qq-y final states contribute if the photon 
is clustered with a quark into a jet and the event is no longer removed by the photon cuts. 

tin Ref. [24] the definition of <5evv was somewhat different and, in particular, did not explicitly contain the effect of 
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2.2 Particle identification 



One of the virtues of e + e~ colliders is the precise knowledge of the energy of the initial state. However, 
ISR of photons can lead to difficulties in the determination of the total energy of the final state. 
Therefore event-selection cuts have been devised to suppress effects due to ISR. In the following we 
describe the procedure employed by the ALEPH collaboration at LEP [53], which we use in our 
numerical evaluations. 

First, particles are clustered into jets according to the Durham algorithm with y cut ,D = 0.002 and 
_E-scheme recombination. Jets where the fraction of energy carried by charged hadrons is less than 10% 
are identified as dominantly electromagnetic and are removed. In the next step the remaining particles 
are clustered into two jets and the visible invariant mass M v - ls of the two-jet system is calculated. 
Using total momentum conservation the reduced CM energy \fp is calculated. The event is rejected 
if s' / s < 0.81. This two-step procedure is later referred to as hard-photon cut procedure (note that it 
is called anti-ISR cut procedure in Ref. [53]). 

Removing events where the photonic energy in a jet is higher than a certain value, as it is done 
in the hard-photon cuts, causes potential problems when perturbatively calculating EW corrections. 
There one relies on the cancellation of infrared (IR) singularities between virtual and real corrections 
when calculating an IR-safe observable. Removing events where a photon is close to a final-state 
charged fermion leads to non-IR-safe observables and spoils this cancellation in the collinear region. 
This feature is common to all observables with identified particles in the final state, and IR-finiteness 
is restored by taking into account a contribution from fragmentation processes, in our case from the 
quark-to-photon fragmentation function. 

3 Structure of the calculation 

To obtain the full EW 0(a 3 a s ) corrections to normalised event-shape distributions and jet cross 
sections we need to derive the NLO EW corrections to the total hadronic cross section and to the 
three-jet production cross section. The total hadronic cross section is decomposed as: 

<7had = / d^ o r™ + / d^L^W + J dafj ^ , (3.1) 

where the first and second terms are the Born and one-loop EW contributions to the process e + e _ — > qq, 
while the last term is the real radiation contribution from the process e + c~ — > qqj. Likewise, we 
decompose the total cross section for three-jet production according to 

J da = J d^ n -^ + / da^™ 

4- [ nV c+c ~^ 9?9 4- [ r |n- e+e ~~*' 9 9 7 4- / r l rr 0+c_ ^9997 , f e + e~ ->qqqq (o o) 

^ / UO virtual,EW ^ / UO virtual,QCD ^ / UO rcal ^ / ""interference ' V°' Z J 

where the first and third terms are the Born and NLO EW contributions of the process e + e~ — > qqg, 
the second and fourth terms are the Born and one-loop QCD contributions of the process e + e~ — > qqj, 
the fifth term is the contribution from the real radiation process e + e~ — > qqg"f, and the sixth term 
results from the contribution of the real radiation process e + e~ — »■ qqqq with identical quark flavours. 
In this work, we are interested in the virtual and real radiation corrections of O (a 3 a s ), which lead to 
the production of three or four jets, when treating photons and hadrons democratically. For the qqqq 
final state, this order corresponds to the interference of the EW amplitude with the QCD amplitude. 
We do not include the squares of the EW and the QCD amplitudes for this process. The former is of 
O (a 4 ) and thus beyond the considered accuracy, the latter is part of the NLO QCD corrections not 
considered in this work. 

We have performed two independent calculations each for the virtual and real corrections, the 
results of which are in mutual numerical agreement. 
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Figure 2: Lowest-order diagrams for e+e — » qq. 




3.1 Conventions and lowest-order cross section 

At the parton level we consider the processes 

e+(fci,o-i) +e~(fc 2 ,(72) -> g(fc 3 ,cr 3 ) +g(fc 4 ,(T4), (3.3) 
e+(fci,CTi) +e _ (fc 2 ,(T2) 9(^3,0-3) + g(fc 4 ,<T 4 ) +g(fc5, A), (3.4) 
e+(fci,o-i) +e _ (fc 2 ,(T2) -> 9(^3,0-3) + g(fc 4 ,CT 4 ) +7(^5, A), (3.5) 

where q can be an up, down, charm, strange, or bottom quark. The momenta k{ of the corresponding 
particles as well as their helicities crj and A are given in parentheses. The helicities of the fermions 
take the values <7j = ±1/2, and the helicity of the gluon or the photon assumes the values A = ±1. 
We neglect the masses of the external fermions wherever possible and keep them only as regulators of 
the mass-singular logarithms. Therefore all amplitudes vanish unless o\ = —o 2 and <r 3 = —0-4, and we 
define a = <j 2 = — G\ and a' = cr 3 = — <j 4 . 

For later use, the following set of kinematical invariants is introduced: 

s = (k 1 +k 2 ) 2 , Sij = (hi + kj) 2 , s ai = s ia = (k a - k { f , a = 1,2, i,j = 3,4,5. (3.6) 

For later convenience we employ the convention that indices a, b = 1, 2 refer to the initial and indices 
i, j = 3, 4, 5 to the final state, while the generic indices I 7 J~ 1, . . . , 5 label all external particles. The 
tree-level Feynman diagrams contributing to the process (3.3) are shown in Fig. 2, those contributing 
to the process (3.4) in Fig. 1, and the ones contributing to the process (3.5) in Fig. 3. 
The lowest-order partonic cross section for the processes given in (3.4) and (3.5) reads 

J da Bo rn - y/c Yl I^ 1 - 2P ^ {1 + 2P2ff ) / d$ 3 \M°/ X \ 2 6 cut ($ 3 ), (3.7) 

A=±l 

where Fq is a colour factor, Pi, 2 are the degrees of beam polarisation of the incoming c + and e~, 
Mq* 7 x is the colour-stripped Born matrix element of the respective process, and the integral over the 
three-particle phase space is defined by 

/ d *3 = (n / jSk?) {27T)4s { ki +k2 -p>)- m 

For the process (3.5), Fq — 3, and for the process (3.4) Fq = 4. The dependence of the cross section on 
the event-selection cuts is reflected by the step function 9 cut ($ 3 ). For the lowest-order cross sections 
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Figure 4: Contributions of all possible vertex functions to e+e — > qqg. 



of (3.4) and (3.5) and the virtual corrections, O cu t depends on three-particle kinematics. It is equal to 
1 if the event passes the cuts and equal to otherwise. 

The formula corresponding to the process (3.3) can be obtained from (3.7) by omitting the depen- 
dence on and the sum over the polarisation A of photon or gluon, using only the two-particle phase 
space $2, and setting Fq = 3. 



3.2 Virtual corrections 

3.2.1 Survey of diagrams and setup of the loop calculation 

We calculate the one-loop EW corrections to the processes given in (3.3) and (3.4), and the one-loop 
QCD corrections to the process given in (3.5). For (3.4) and (3.5) their contributions to the cross 
section are generically given by 



/ 



da virtU ai = y/c j(l-2 J PH(l + 2P 2( T)y"d$3 2Re[Mr' A (^r' A )*]ecut($3), (3.9) 



a,cr — ± 
A=±l 



where the notation is the same as in the previous section and Ai1 a x denotes the contributions of 
the virtual corrections to the matrix element after splitting off the colour factor of the corresponding 
lowest-order amplitude. 

The NLO EW virtual corrections to (3.3) and (3.4) receive contributions from self-energy, vertex, 
box, and in the case with a gluon in the final state, also pentagon diagrams. The structural diagrams 
for the process with gluon emission containing the generic contributions of all possible vertex functions 
are shown in Fig. 4. The structural diagrams for the process without gluon emission can be obtained 
by taking the first four and the sixth diagrams of Fig. 4 and discarding the outgoing gluon. The 
Feynman diagrams contributing to the 3-point vertex functions are shown in Fig. 5, those contributing 
to the 4-point and 5-point vertex functions in Fig. 6. The diagrams for the Z-boson, photon, and quark 
self-energies can be found for example in Ref. [54]. The symbol q stands for the quarks appearing in 
(3.3)-(3.5), the symbols q' for their weak-isospin partners. Since we neglect the masses of the external 
fermions wherever possible, there are no contributions involving the physical Higgs boson coupling to 
those particles. For b quarks in the final state, diagrams with W bosons also have counterparts where 
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the W bosons are replaced by would-be Goldstone bosons, which are included in the calculation but 
not shown explicitly in the figures. We also do not depict diagrams that can be obtained by reversing 
the charge flow of the external quark lines in the first six diagrams of Fig. 4. 

In total we have 0(200) contributing diagrams in the 't Hooft-Feynman gauge for the process 
with gluon emission and O(80) for the process without gluon emission, counting closed- fermion- loop 
diagrams for each family only once. 

The NLO QCD virtual corrections to (3.5) receive contributions from self-energy, vertex, and box 
diagrams. The corresponding Feynman diagrams are shown in Fig. 7, where we have omitted quark 
self-energy contributions. We do not depict diagrams that can be obtained by either reversing the 
charge flow of the external lepton lines in the first diagram or of the external quark lines in the last 
three diagrams of Fig. 7. In total we have O(20) contributing diagrams in this case. 




Figure 7: Sample diagrams for virtual QCD corrections to the process e+e — ► qq-f. 
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We treat the gauge-boson widths using the complex-mass scheme, which has been worked out at 
the Born level in Ref. [55] and at the one- loop level in Ref. [56]. In this framework the masses of 
the Z and the W boson are complex quantities, defined at the pole of the corresponding propagator 
in the complex plane. As a consequence, derived quantities like the weak mixing angle also become 
complex, and the renormalisation procedure has to be slightly modified. Introducing complex masses 
everywhere in the Feynman rules preserves all algebraic relations like Ward identities and therefore 
also gauge invariance. Terms that break unitarity are beyond the one- loop level (see Ref. [56]). 

3.2.2 Algebraic reduction of spinor chains 

We have performed two independent calculations of the virtual corrections. In version 1 of our cal- 
culation we generated the amplitudes using FeynArts 3.2 [57] and employed FormCalc 5 [58] to 
algebraically manipulate the amplitudes, which led to 150 different spinor structures. In order to re- 
duce the number of spinor structures, we applied the algorithm described in Ref. [56] and extended 
it to the case with one external gauge boson. In this way, we reduced all occurring spinor chains to 
O(20) standard structures, the standard matrix elements (SMEs), without creating coefficients that 
lead to numerical problems. 

After the reduction of the spinor structures, we separate the matrix elements into invariant co- 
efficients F n and SMEs M. The F n are linear combinations of one- loop integrals with coefficients 
depending on scalar kinematical variables, particle masses, and coupling factors, the M contain all 
spinorial objects and the dependence on the helicities of the external particles (see e.g. Ref. [59]): 

M oo'\ = J- F ™' x {{s, Sij ,s ai }) M a / X (fci, A&, fo, fed, h) . (3.10) 

n 

All contributions to the matrix element involve the product of two spinor chains corresponding to the 
incoming leptonic and the outgoing hadronic current. These spinor chains can be contracted with one 
another, with external momenta, or with the polarisation vector of the outgoing photon or gluon. 

In the following we describe the strategy to reduce all occurring products of spinor chains and 
polarisation vectors to a few standard structures. In this section we choose all particles incoming and 
use the short-hand notation 

[A]± b = Va(ka)Au±u b (k b ) (3.11) 

for a spinor chain, where v a {k a ) and ut(fet) are spinors for antifermions and fermions, respectively, 
with the chirality projectors to± = (1 ± 7s)/2. We denote the external polarisation vector by e. Since 
we work with massless external fermions, only odd numbers of Dirac matrices occur inside the spinor 
chains. 

The objects we want to simplify are of the form 

Au p u 2 {k 2 ) x v 3 {k 3 )Buj r u 4 {k 4 ) = [A] P 12 [B] T 34 . (3.12) 

We make use of the Dirac algebra, the Dirac equation for the external fermions, transversality of the 
polarisation vector, momentum conservation, and relations resulting from the four-dimensionality of 
space-time, which can be exploited after the cancellation of UV divergences, which are dimensionally 
regularised in our work.* In four dimensions one can relate a product of three Dirac matrices to a sum 
where each term only consists of a single Dirac matrix multiplied by either the metric tensor or 75 
and the totally antisymmetric tensor t lxvpa through the Chisholm identity as 

VW = g^Y - a w i v + g" p r + ^ vpa i^- (3.13) 

A further consequence of the four-dimensionality of space-time is the fact that one can decompose 
the metric tensor g^ v in terms of four linearly independent orthonormal basis vectors n; (see e.g. 

■tlR divergences, which may be regularised dimensionally, do not pose problems in this context, as shown explicitly in 
the appendix of Ref. [60] . This fact is confirmed in our calculations where we alternatively used a four-dimensional mass 
regularisation scheme or dimensional regularisation, leading to the same IR-finite sum of virtual and real corrections. 
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Refs. [61,62]) 



(3.14) 



k,l=l 



1=1 



with (nkni) = gki, where (gki) = diag(l, — 1, — 1, — 1). A convenient choice of the four vectors n; in 
terms of three linearly independent massless external momenta fcj, kj, kk is given by 



n (ki, fcj) fcfc) = 

Tig (fej, fcj, fcfc) = 



1 



v / 2(k l k :j )(k l k k )(k : jk k ) 
1 

v /2(k i k J )(kik k )(k : jkk) 



(kjkk)k^ -\- (kikk)kj' (kikj)k^ 



T h h- hi 



(3.15) 



In particular, the construction of the forth independent momentum via the totally antisymmetric 
tensor e^'* 7 avoids the appearance of inverse Gram determinants. 

For the reduction of spinor chains we use the following algorithm. In the first step we disconnect 
two spinor chains which are contracted with each other using the decomposition (3.14), 



la ® 7 M = 



i - • 



(3.16) 



i=i 



The choice of the external momenta in the above decomposition strongly depends on the position of 
the contracted Dirac matrices inside the spinor chain. It is advantageous to choose them in such a way 
that one can make use of the Dirac equations v(ki)fii = and h'iuiki) = and the mass-shell condition 
fif = kf = in a very direct manner, avoiding unnecessary anticommutations. We follow the algorithm 
described in detail in Ref. [56]. After simplifying the expressions using the identities above, there are 
remaining contributions from the contraction of a basis vector n 3 with a Dirac matrix inside the spinor 
chains, which can be eliminated by employing the Chisholm identity as 



'^■{fiifijfik {kikj)$k ~t~ {kikk)$j [kjkk)^i\ 75 

y/2{kikj){kikk){kjkk) 



(3.17) 



In the calculation at hand, we have to deal with a maximum of three contractions between the two 
spinor chains. After disconnecting them, we are left with objects of the form [f\^ b , where the vector 
p can either be an external momentum kj, j ^ a, b, or the polarisation vector e of the external gluon 
or photon. 

In the next step, we reduce the spinor chains that do not involve e using the relation 



(3.18) 



i=i 



Choosing the indices i, j, k in (3.15) appropriately, this allows to eliminate all external momenta in the 
spinor chains apart from one for each chain via 



fi'n 



(3.18) 



ab 



km.fi ^ 1 9 n j{^ai fcb? fcn) ^i(&aj &n) 



(3.17) (kgk n ) (k b k m ) - (kgk b ) (k n k m ) + (kgk m ) (fc n fc t ) ± i ^ " 0° kg , g fc„ ,„ fcfc , g fc m , g 

2 (kgk n ) (kbkn) 



fin , (3.19) 
- ab 



where m ^ a,b, n. The described reduction allows us to express all occurring spinor structures in terms 
of a linear combination of 20 SMEs 



e]? 2 [*i]»4» N12 H34 5 (e fc i)[ fc 3]i 2 [A:i]34, (ek 2 ) [k 3 } a 12 [ki} r 34 , (ek 3 ) [k 3 } a 12 [ki} T 34 . 



(3.20) 
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The reduction to this basis introduces at most two summands per spinor chain. Inserting the SMEs 
(3.20) into the amplitude reduces its size by a factor of two. Since different reduction strategies did 
not lead to more concise results we chose to use the SMEs given in (3.20) in this calculation. 
For the virtual corrections to Chad, where k$ and e are absent, we use the four SMEs 

MI2 N34 ■ (3-21) 

Version 2 of our algebraic calculation starts from diagrammatic expressions for the one-loop cor- 
rections generated by FeynArts 1.0 [63] and proceeds with the algebraic evaluation of the loop 
amplitudes with an in-house Mathematica program. The algebraic manipulations do not make use 
of four-dimensional identities for the Dirac chains for the 2^3 process, so that a larger set of 64 
SMEs had to be introduced. After rendering the corresponding invariants F n UV finite upon adding 
the counterterms from the renormalisation and IR finite upon including the "endpoint contributions" 
from the subtraction function (see Section 4.2), the SMEs are evaluated in four space-time dimensions 
using the spinor formalism described in Ref. [64]. For the 2 — > 2 process the only Dirac structure that 
involves divergences is [7 M ]i2bVi]34> and all other SMEs are proportional to these, as can be deduced 
from the identities given in (3.10) of Ref. [56]. 

3.2.3 Evaluation of the loop integrals 

The tensor integrals are evaluated as in the calculation of Refs. [56,65], i.e. they are numerically reduced 
to scalar master integrals. The master integrals are computed using complex masses according to 
Refs. [66-68], using two independent Fortran implementations which are in mutual agreement. Results 
for different regularisation schemes are translated into each other with the method of Ref. [69] . Tensor 
and scalar 5-point functions are directly expressed in terms of 4-point integrals [70,71]. Tensor 4- 
point and 3-point integrals are reduced to scalar integrals with the Passarino-Veltman algorithm [72] . 
Although we already find sufficient numerical stability with this procedure, we apply the dedicated 
expansion methods of Ref. [71] in exceptional phase-space regions where small Gram determinants 
appear. 

UV divergences are regularised dimensionally. For the IR, i.e. soft or collinear, divergences we either 
use pure dimensional regularisation with massless gluons, photons, and fermions (except for the top 
quark), or pure mass regularisation with infinitesimal photon, gluon, and small fermion masses, which 
are only kept in the mass-singular logarithms. When using dimensional regularisation, the rational 
terms of IR origin are treated as described in Appendix A of Ref. [60] . 

3.3 Real Corrections 

In this section we describe how we evaluate the last two terms in (3.2). The real radiation corrections 
to the total hadronic cross section [last term in (3.1)] are computed along the same lines. 
The processes we consider are given by 

c + (fci, <7i) + e~(/c2,er 2 ) -> q(k 3 ,a 3 ) + g(/c 4 ,er 4 ) + g(fc 5 , Ai) + "f(k 6l A 2 ), q = u, d,c,s,b, (3.22) 
c + (fci, (Ti) + e~(/c 2 , 02) q(k 3 ,a 3 ) + q{ki, 1T4) + q(k 5l a 5 ) + q(k 6 ,a 6 ), q = u, d, c, s, b. (3.23) 

The corresponding cross section is obtained as 

fda ieal = ± F c ]T \(l-2P 1 a)(l + 2P 2 a)\ ]T f d$ 4 \M°£%*\ 2 6 cut ($ 4 ) 

+ £ /d$ 4 2Rc{(^;^)^^ CD }e cut ($4)), (3.24) 

where we have used helicity conservation to simplify the helicity sum, such that a denotes the helicity 
of the incoming electron, and a' and a" = <r 5 = — er 6 the helicities of the outgoing quarks. The 
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Figure 8: Sample diagram for a non-trivial interference between EW and QCD diagrams in e+e — > qqqq. 



four-particle phase-space volume element reads 

1 d 3 fc3 d 3 fc 4 d 3 fc 5 d 3 fc 6 . , 4 „ , , , , , . . n . 

4 = (2^ w W W W ( } ( 2 ~ " 4 " ~ } ' ( } 

and i<c = 4 for both contributions. As in (3.7), C ut ($4) represents cuts used in the event selection. 

We only consider the 0(a 3 a s ) corrections in our analysis. At this order, the process (3.23) receives 
a contribution only from the interference of EW with QCD amplitudes, as illustrated in Fig. 8. Owing 
to colour conservation, this interference term is only non-zero for identical quark flavours and for 
a 1 = a" . It is non-singular over the full phase space defined by the event-selection cuts. 

The integral of the process (3.22) over the four-particle phase space contains IR divergences due 
to the emission of a soft or collinear photon or gluon. Prior to numerical implementation, one has 
to isolate these divergences and combine them with the corresponding contributions from the virtual 
corrections. In our implementation, we use three different methods for this task: two variants of 
the phase-space slicing method [73-79], and the dipole subtraction method [80,85,86]. In phase- 
space slicing, the phase space is split into singular and non-singular regions. In the singular regions 
the integration over the singular variables is performed analytically, while the non-singular regions 
are integrated over fully numerically. In the dipole subtraction method, a subtraction function that 
mimics the singular behaviour of the integrand is added and subtracted, leaving a finite four-particle 
phase-space integration and a remainder, where the integration that leads to singularities is carried 
out analytically. Both methods rely on factorisation properties of the matrix elements and of the phase 
space in the soft and collinear limits. The singularities are treated analytically so that a numerical 
integration does not pose any problems. 

Both methods described above are valid for NLO calculations. We employ them separately for the 
photon and the gluon in the calculation at hand. However, they are not sufficient in the region where 
both the photon and the gluon become soft or collinear at the same time. This region corresponds to 
two-jet production where the photon and the gluon are unresolved. In this region close to the kinematic 
endpoint in the event-shape distributions fixed-order calculations are in any case not appropriate. To 
prevent problems with numerical stability from this region, we impose a lower cut-off on each event- 
shape variable in the first bin of the distribution that corresponds to two-jet production. 

To be able to compare the results of our calculation to experimental measurements and to improve 
the accuracy of the theoretical prediction relevant for the a s determination, we have to incorporate 
the kinematical cuts used in the specific experiment. For the event-selection cuts, C ut, we apply 
the procedure used by the ALEPH experiment, as described in Section 2.2. Electromagnetic jets are 
removed by imposing an upper cut on the photon energy in the jet. However, the cut on the photon 
energy also removes events with a highly energetic photon collinear to a quark or antiquark in the 
final state that lead to a configuration where the photon and quark are in the same jet. The left- 
over collinear singularity associated with this isolated-photon rejection is properly accounted for by a 
contribution from the quark-to-photon fragmentation function (see next section). 

Collinear photon emission in the initial state is regulated by the mass of the electron and leads to 
large contributions of the form a n ln n (m^/s). The cut on the visible invariant mass M v - ls removes part 
of the hard photon emission collinear to the incoming beam particles and therefore suppresses large 
mass-singular logarithms. 
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4 Treatment of soft and/or collinear singularities 

In this section we describe the treatment of the soft and collinear singularities for the process (3.22). 
4.1 Phase-space slicing 

In the phase-space slicing method, technical cut parameters are introduced to decompose the real 
radiation phase space into regions corresponding to soft or collinear configurations (unresolved regions), 
and a resolved region that is free from singularities. Consequently, the cross section decomposes into a 
soft, a collinear, and a finite part. Since we apply an isolation cut on the final-state photons (i.e. on a 
specific, identified final-state particle), we must also include a fragmentation contribution, so that the 
total real contribution reads 



d(T rcal = d<7 so f t + dcr C oll — d(7f rag + definite- (4.1) 

In the following sections we first review the method for collinear-safe observables, the necessary mod- 
ifications for non-collinear-safe event-selection cuts are described in Sections 4.1.2 and 4.3. 



4.1.1 Collinear-safe observables 

Soft and collinear contributions contain IR divergences which are evaluated analytically, while the 
finite contribution is evaluated numerically. Since no quark-mass singularities from final-state radi- 
ation remain for collinear-safe observables, no fragmentation contribution is necessary in this case, 
i.e. do-fr a g = in (4.1). We implemented two different variants of phase-space slicing: the two-cut- 
off slicing according to Refs. [79,81], which uses mass regularisation, and the one-cut-off slicing of 
Refs. [76,77] within dimensional regularisation. The application of both methods is described in detail 
below for IR-divergent contributions due to unresolved photons. The unresolved gluon contributions 
are obtained by an appropriate replacement of the coupling constants. 

In two-cut-off slicing, the splitting of the phase space into singular and non-singular parts is achieved 
by introducing a cut S s on the photon energy E 1 < d s y/s/2 — AE in the CM frame. The collinear 
region is defined by E 1 > AE and 1 > cos 8 > 1 — S c , where 6 is the smallest angle between the photon 
and any charged fermion in the CM system. 

In the soft region the squared matrix element |A4 roa i| 2 and the phase-space element d$4 factorise 
so that we can apply the soft-photon approximation, e.g. given in Refs. [59,82], and introduce an 
infinitesimal photon mass ra 1 as a regulator. The resulting soft-photon contribution is 



da soft = da Born f V V (-l) I+J Q I Q.j (2 In (—) [2 - In ( 
Z7T 1=1 j=j+i l V m 7 / L 



2 In 



UEjEA 


n ln2 


f4Ej\ 




\ mimj J 









m . 



Hze! 
3 



2 Li 2 1 



AE I E,j 
sij 



(4.2) 



where we keep the masses of the fermions as regulators for the collinear singularities (Ei 3> m/), i.e. we 
have mi = m 2 = m c and m 3 = TO4 = m q , and dcrBorn denotes the Born cross section for e + e~ — > qqg. 

Hard collinear contributions arise from the limit where the photon becomes collinear either with one 
of the incoming beams or with the outgoing (anti-)quark. These contributions contain mass-singular 
logarithms, which are regularised by the fermion masses. Their integrated forms are again proportional 
to the Born cross section for e + e~ — > qqg. We split the collinear photon contributions into those from 
initial-state radiation and those from final-state radiation: 



dcr, 



coll 



j initial , j final 
acr coll + a<7 coll ' 



(4.3) 



In the case of initial-state photon emission, the available e+c~ CM energy is reduced, and the Born 
process is probed at this reduced energy. Moreover, the mass regularisation introduces a spin-flip term. 
We indicate these two features in the argument of the Born cross section, where ki and Pi being are 
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the momenta and the degrees of polarisation of the respective beams. The integrated initial-state and 
final-state collinear contributions thus read 



j initial 
dcr coll 



= E 



0=1 



27T^ a 



1-6. 



<{ dfTBorn (ka ~> zk a ) Pff{z) In _ 1 

+ dfTBorn (fca 2fc a , P a -> -P a ) (1 - z) 



and 



where 



4 , 

do-col? 1 = ^Tj^Qi dcr Born | 



2 +2 "{-W) 



1 - In 



2£ 2 <S C 



3 - 



2tt 2 



1 + 



(4.4) 



(4.5) 



(4.6) 



denotes the / — > /7 splitting function. The parameters <5 S and <5 C govern the splitting of the phase space 
in the different regions, but the final result does not depend on them. They have to be chosen small 
enough to guarantee that the applied approximations are valid. Therefore, varying these parameters in 
a certain range and showing independence of the results serves as an important check of the calculation. 

In the one-cut-off slicing, soft and collinear regions are defined by a cut y m ; n on two-particle 
invariants su = (kj + kj) 2 normalised to the e + e~ CM energy squared, yjj — Su/s. The soft region 
of parton / is defined by \yu\ < y m in for all J, while the collinear region of partons / and J is defined 
by \yu\ < 2/min, while all other invariants \yix\ > J/min for K ^ J. 

Following Ref . [77] , the soft and collinear divergent contributions are first calculated in the unphys- 
ical kinematical situation of all particles outgoing, and later on continued to the physical kinematics 
by means of a crossing function, such that 



do-coil 



A ~ out i A ^cross 
d(7 coll + dcr o 



'coll 



(4.7) 



The combination of the first term with the soft contribution yields [76] 

4 

d^oft+d^ = daBom^j^^y £ (-1) I+J QiQj 



I, J=l 



\su\ 



3 



\su\ 



In 2 



\yij\ 

ymin 



3, 
+ 2 ln 



\yij\ 



7 7T 



2 1 



(4.8) 



where we used dimensional regularisation in d = 4 — 2e dimensions with mass parameter \x required to 
maintain a dimensionless coupling. 

Kinematical crossing of c+c~ to the initial state introduces crossing functions, which were derived 
for QCD using factorisation of parton distributions in the MS scheme in Ref. [77]. For photon radiation 
off incoming electrons, this MS expression is converted using Ref. [79] to a mass-regularised form, 
involving the physical electron mass. This results in the following crossing term, 



j cross 
Q(T coll 



V 2tt 

a=l,2 



Ql / dz do-Born {k a 
JO 




S(l-z) 



-1 \ +l-z 



+ [p // (*) (ln (ymin)+ ln(-^— ^ 

where we use the usual [. . .] + prescription 

/ dx [f{x)} + g{x) = [ dxf(x) [g{x) - g{l)} . 
Jo Jo 



(4.9) 



(4.10) 
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The results (4.8) and (4.9) apply to unpolarized cross sections, but the polarization effects of the 
incoming particles can be easily restored by the spin-flip terms of (4.4), 

pi 

dcr co°ll = dcr coll + I dz d<J Born {k a -> zk a , P a -> - P a ) ~ der B orn (k a -> zk a ) (1 - z) . (4.11) 



o=l,2 



4.1.2 Non-collinear-safe observables 

The upper cut on the photon energy inside jets affects the slicing procedure. Imposing this cut in 
the soft and finite parts of (4.1) is trivial. Since the cut acts outside the soft region, this part is not 
affected at all. Only the collinear singular part needs some non-trivial adjustments. More precisely, only 
collinear final-state radiation needs to be considered, since photons collinear to the incoming electrons 
or positrons can never appear inside jets within detector coverage. In Section 4.1.1 we parametrised 
the collinear region in terms of the energy fraction z carried by the (anti-)quark that results from the 
splitting. The experimental cut, however, is a cut on the energy fraction 1 — z of the photon after the 
splitting, i.e. the cut on the energy fraction of the photon reads 1 — z < z cut . 
In the two-cut-off slicing approach, imposing this cut generalises (4.5) to 



4 /.l-AE/Ej C 

da^\z cut ) = ^-Q>da Boin / dz P ff (z) 



i=3 



In 



1 



+ (l-z) . (4.12) 



Performing the integration yields 

4 

dCT colf (^cut) = X! dCTBorn ' 



i=3 



-2z cu t + 



'cut 



2 In 



fAE/Ei 



Zcut 



In 



f2EfS c 



\ 2 



2 In 



fAE/EA 



- 4Li 2 (zcut) + (1 - zcut) (3 - zcut) ln(l - z cut ) + 5z cut - 



Zcut / 

(4.13) 



In (4.13) we have the original dependence on the slicing parameters and the mass regulators, already 
present in (4.5), plus an additional term that depends on the slicing parameters, the mass regulators, 
and the cut on the photon energy. It is exactly this term that gives rise to left-over singularities. 

Within the one-cut-off approach, a cut on the final-state photon energy in jets can be conveniently 
taken into account by including a collinear photon contribution in d<7 co ii, 



do-coil = da-^* + d(T^°j ss - dcrJ oll (z cut ), 



(4.14) 



which subtracts the contributions of collinear photons with energies above z cut . This contribution 
reads 



dcrJ oll (z cut ) = 



4 

2^ 

i=3 



Qid(7Born 



i=3 



— QidfTBorn 




where the e-dependent splitting function pjy* is given by 



Pf(z) 



1 - z 



(4.15) 



(4.16) 



The derivation of this collinear unresolved photon factor is described in detail in Refs. [83,84]. 

In both slicing approaches, we observe that the introduction of the cut on the final-state pho- 
ton energy results in uncompensated collinear singularities in (4.13) and (4.15). These are properly 
accounted for by factorisation of the quark-to-photon fragmentation function, explained in detail in 
Section 4.3 below. 
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4.2 Dipole subtraction 



The basic idea of any subtraction method is to subtract an auxiliary function from the integrand that 
features the same singular behaviour in the soft and collinear limits. The partially integrated auxiliary 
function is then added to the virtual corrections (and counterterms from the factorisation of parton 
distributions and fragmentation functions) resulting in analytic cancellation of IR singularities. We 
use the dipole subtraction method, first introduced in Ref. [80] for massless QCD and later generalised 
to massive fermions for collinear-safe and non-collinear-safe observables in Ref. [85] and Ref. [86], 
respectively. Since we regulate IR divergences with particle masses, we follow the description of 
Refs. [81,85,86]. Again we first describe the method for the case of an unresolved photon; the gluon 
case is obtained by a suitable replacement of couplings. In this section we suppress the weighted sum 
over initial-state polarisations. 



4.2.1 Collinear-safe observables 



Suppressing flux and colour factors, the 6-function related to the phase-space cuts, dipole subtraction 
is based on the general formula 



J d$ 4 ^ l-Mrcall 2 = J d$4 ^ |Xrcal| 2 - \M suh \ 



d$ 4 |X S ub| 2 , 



(4.17) 



where A labels the photon polarisation. The subtraction function M S ub is constructed from ordered 
pairs IJ of charged fermions, where fermion / is called the emitter and fermion J the spectator. Only 
the kinematics of the emitter leads to singularities. The spectator fermion is used to balance energy- 
momentum conservation when combining the momentum of the photon with the momentum of the 
emitter. Making the dependence of the subtraction function on the emitter-spectator pair explicit, we 
can write 



\M 



IJ, sub 



($ 4 )| 



-(-1) 



I+J 



,/J.Tcr/J 

T=± 



(4.18) 



where aj is the helicity of the emitter and k = k e the photon momentum. The sum over r = ± takes 
care of a possible spin-flip resulting from collinear photon emission. The subtraction function A4i,j, su b 
depends on the whole four- particle phase space $4, whereas the Born matrix elements depend only 
on three-particle phase spaces $3,77. The mappings of the four-particle on the three-particle phase 
spaces, which are different for each emitter-spectator pair IJ and explicitly given in Ref. [85], ensure 
proper factorisation in each singular limit. 

In the massless case, the dipole factors are explicitly given by 

9ij i+ \h, kj,k) 
9ia i+ \ki, k a , k) 
9 a i,+ {ka, ki, k) 
9ab,+ (ka, kb, k) 



(hk) (1 - yij ) 
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{k a k)x a b 



2 Xia %it 
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2-xi 
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/-1 \ 1 Z ij 



1 %io 
- 1 - Xi 



1 - x ab 



-l-Xab 



(sub) , 



g>-;-'(ki, k jt k) = gt^ih, k a , k) = g { :£\k a ,k u k) = g a T*\k a , k b , k) = 0, (4.19) 

where we denote final-state particles with the letters i, j and initial-state particles with the letters a, b 
and use the definitions 

kdkfy kftk ■ k\)k k^ki ~\~ k^k ■ k%k kjk 

%ab = ; i ; %ia i j ; J j ; Uij 



ka kb 



kaki -\- k a k 



ki kj 



ki k ~\~ kjk 



ka ki 



ki kj 



kaki k a k 



ki k j ~\~ kjk 



(4.20) 
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The finite part of the real corrections thus reads 



d$ 4 



©cut($3,/j) 



I.. 7=1 



(4.21) 



where from now on we suppress the spin sums of the fermions in the notation. The cuts on the 
momenta of the final-state particles are included in terms of the functions O cu t(^4) and @ C ut(^3,-Z"j)> 
where the arguments signal which momenta are subject to the cuts. Note that already at this point 
collinear safety is assumed, because the emitter momentum in $3,77 tends to the sum of emitter and 
photon momenta in the collinear limit by construction, i.e. a recombination of these two particles in the 
collinear limit is understood. This has to be changed in the treatment of non-collinear-safe observables 
considered below. Note that the vanishing of the spin-flip parts gi.j.- only holds in the difference (4.21), 
where the collinear singularities cancel; the spin-flip parts contribute when the collinear singular region 
in l-M/j^ubl 2 is integrated over, as will become apparent in the following. 

We turn now to the treatment of the singular contributions. Splitting the four-particle phase space 
into a three-particle phase space and the photonic part, 



J d$ 4 = J dxj d$ 3 ,/j(aO J 



d$ 



1JJ1 



(4.22) 



where the photonic part of the phase space, d$ 7i /j, depends on the mass regulators mj and m 7 of 
the fermions and of the photon, we have to compensate for the reduction of the CM energy in the 
case of ISR radiation, which is indicated by the x-dependence of the three-particle phase space. The 
integration over x, thus, becomes process dependent and cannot be carried out analytically in general. 
Instead, it is possible to split off the singular contribution resulting from the endpoint of this integration 
at x — > 1 upon introducing a [...]+ distribution. Leaving details of the analytic integrations over the 
singular phase spaces d$ 3 jj(x) to Ref. [85], the result for the integrated singular part of the real 
corrections reads 



/ 



Z7T i,J=i r =± 



+ Ysj d<I '3GS ) ( S /./)|XBorn($3,T ( 7 / )| 2 e cut ($3)|, 



(4.23) 



where the functions Gj S j^ (sjj, x) and G^ff^) (s/j) result from the analytic integration over the photonic 
part of the phase space. For the final-final emitter-spectator case, there is no convolution part Q }| , 



i.e. 



g£V(a ii ,x)=0, G^\ Sij ) =8^jd^ ij g^(k i ,k j ,k), (4.24) 



in contrast to the other emitter-spectator cases, where 



G^\sij,x) = 87t 2 x J d^jjgfffikj,^,^, G^( SlJ )=j\xg^ 



( SlJ ,x). (4.25) 



consistently takes 



It should be noted that the invariant §u in the integration over 

the values Sjj = Zkjkj, i.e. in the "endpoint" at x = 1 this variable corresponds to the invariant 
sjj = 2kjkj of the three-particle phase space <J> 3 = $3jj(x = 1) corresponding to the phase space 
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without photon. The explicit results for the functions Q and G read 



and 
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l-x 
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QSHS ai ,x) =Pff(x) 
Ga£+(§ ab ,x) = Pff(x) 



■ln(2-a;) + (l+a;)ln(l-a;), 



gfSis^x) = gi™°>(~s ia ,x) = o, g';c>(s ab ,x) = g^_\s ai , x ) = i - 



- rS suh hz 



(sub) ,~ 



- rM uh hz 



= i - x 



^ + 1, G^\s ia )=C(\s ia \,' l 



G$?(8ai) = £ (\sai\,m 2 a ) + y - ^ GabX } (sab) = C (s ab ,m 2 a ) + 
G^\ Sij ) = G^(s ia ) = G^( Sai ) = G ( :t_\ Sab ) = 1 
where the auxiliary function 



7T 2 * 



In ^ 



Inl^i 



(4.26) 



(4.27) 



(4.28) 



contains the soft and collinear singularities of the endpoint parts. These G terms, which have the same 
kinematics as the lowest-order contribution, exactly cancel the corresponding singular contribution 
from the virtual corrections. 



4.2.2 Non-collinear-safe observables 

In Ref. [86] the dipole subtraction method for collinear-safe photon radiation [85], as briefly described 
above, has been generalised to non-collinear-safe radiation. There, in the collinear photon radiation 
cone around a charged final-state particle, the fraction z of the charged particle's energy is not fully 
integrated over, but the cut procedure affects the range of the z integration. In the auxiliary three- 
particle phase spaces $3^(2;) introduced in Section 4.2.1 the role of z is played by the variables Zij 
and Zi a defined in (4.20). The transition from the collinear-safe to the non-collinear-safe case requires 
both a modification in the subtraction procedure and a non-trivial change in the analytical integration 
of the subtracted parts that will be re-added again. We briefly describe these generalisations in the 
following and refer to the original formulation in Ref. [86] for more details. 

In the subtraction procedure, as given in (4.21), the cut prescription C ut( ( I > 3. ij) is modified in 
such a way that the auxiliary momentum ku for the emitter-photon system is replaced by the two 
collinear momenta k t — Zijku and k = (1 — zu)ku for the emitter and the photon, respectively. This 
procedure concerns only contributions of final-state emitters. 

On the side of the re-added subtraction function this implies that the cut on Zij has to be respected 
as well. In practice this means that in most cases at least the non-singular parts of this integration 
have to be performed illy The detailed procedure is described in the following separately for 

the two cases of final-final and final-initial emitter-spectator pairs. 

Final-state emitter and final-state spectator 

The integration of the radiator functions for a final-state emitter i and a final-state spectator j is of 
the form 

G fc b) (««) = Y J dl/«(l-y«) J d Zl3 g^\h,k 3 ,k), (4.29) 
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where yij and Zij are given in (4.20). The limits of integration, which can be explicitly found in 
Ref. [85], depend on the regulator masses and Lorentz invariants of the emitter, spectator, and photon 
system. The explicit integration leads to the results given in (4.27) in the massless limit. In the non- 
collinear-safe case we want to use information on the photon momentum in the collinear cone, which 
is controlled by the variable z^, i.e. we have to interchange the order of the integrations in (4.29) and 
leave the integration over open. To this end, we consider 



(4.30) 



where the limits 2/1,2(2^) of the depend on the mass regulators (see Ref. [86] for details). The soft 
singularity contained in (4.30) can be split off by employing a [...]+ distribution, 



§^\~ Sij ,z) = G^\ Sij )6(l -z)+ teW*) 



(4.31) 



so that this singularity appears only in the quantity G\^\sij) already given in (4.27). In the limit of 
small fermion masses, the integral in (4.30) can be carried out explicitly, resulting in 



G$i ) (S ij ,z) = P ff (z) 



KSf) 


- 1 







+ (l + z)ln(l-z), 



Gii"-(Sij,z) = 1 - z. 



(4.32) 



The explicit form of the ij contribution \A4 S ub,ij ($i) | 2 then reads 

x [aft? (i« 



)s(i-z)+ g\tr ] {hj,z) 



Mo{ki,kj\T(ii^ Q cut (ki = zki,k = (l-z)ki,kj,{k n }^j . (4.33) 

The term in curly brackets in (4.33) consists of a term proportional to a (^-function in z, which is 
the usual endpoint contribution, and a [. . .] + prescription, acting only on 6 cut . In our case 6 cut just 
provides a lower cut-off on the z-integration, 9 cut = 9{z — 1 + z cut ), and we find 



J d$i|M Bub ,«($i;ai)| 2 - -T^Y,{-l) i+i QiQj f d *o 



Mo ( h, kj\ to, 



x {G\ 



ij.T ^ / 
JO 
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(4.34) 



The z-integration in the second term of (4.34) can be carried out explicitly, and the sum over r can 
be performed, because we consider only unpolarized final states. In this way we obtain 

...1. ■■. i I / s ■ \ 

ln(zcut) 



Si 



IT 

~3~ 



- - 2 In 
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^(1 - Zcut) 
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&ij (1 ^cut) 

m?z cut 



2Li 2 (z C ut) ■ 

(4.35) 



Final-state emitter and initial-state spectator 

In the case of a final-state emitter i and an initial-state spectator a, the integration of |.Mia,sub| 2 over 
x = x ia is performed using a [. . .], prescription, 



J o dx ia J dz ia g { ^\k u k a ,k)--- = dx^G ( - s ^\s ia )S(l-x) + £?j™ r b) (s ia ,x) 



(4.36) 
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where the ellipses stand for x-dependent, process-specific functions like the Born matrix element 
squared or flux factors. Here we used the fact that the upper limit x\ of the x-integration is equal to 
one up to mass terms that are only relevant for the regularisation of the singularities which eventually 
appear in G . The integration over x is usually done numerically. Details of the derivation of the 

functions G[™^ and G^^\ which are given in (4.27) and (4.26), can be found in Ref. [85]. 

In the non-collinear-safe case (see again Ref. [86] for details), the ellipses also involve terms like the 
cut function that depend on Zi a . Therefore the whole integration has to be done in general numerically. 
To this end, a procedure is employed that isolates the occurring singularities in the endpoint. The basic 
idea in this procedure is the use of a generalised [. . .] + prescription that acts on multiple variables. 
Denoting the usual [. . .] + prescription in an n-dimensional integral over the variables r^, i = 1, . . . , n, 

by 

J d"r [ff(r)]g o) /« ee J d"r 5 (r) (/(r) - /(r)| n=a ) , 
the natural generalisation to two variables reads 



(4.37) 



/ d n r[ 5 (r)]^: : »/(r) ee J d»r [fo(r)]<# o) 



(r 3 ) 



= J d"r 5 (r)(/(r)-/(r)| r4=0 -/(r)| rj=6 + /(r)|r J =.). (4.38) 

To recover the usual notation, we drop the subscripts a and/or b if they are equal to one. The generic 
form of the integral we want to perform is 



I[f] = 



22 (x) 

zi(x) 



dzg^\ki,k a ,k)f(x,z), 



(4.39) 



with f(x,z) denoting an integrable function of x = xi a and z = Zi a . After introducing the multiple 
[...] + distribution, all soft and collinear singularities are integrated out, and the result contains only 
regular integrations over x and z within unit boundaries, 



/[/] = f dx f dz 
Jo Jo 
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(4.41) 



Using this result, the explicit form of the ia contribution \M su b,ia ($i) | 2 reads 

J d$i|M Bub ,i„($i;o-i)| 2 = -^ : 2(-l) 0+i Q Q i jr dx J d$ 0iia (s ia ,x) jf dz 

-,(sub) /_ 
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Inserting the explicit form of the cut function, 6 cu t = 0(z — 1 + 2 C ut), this yields 
J d*i \M suh , ia ($i;o-i) | 2 = ~ ]T (- l ) a+l Q*Q* 



d$ , M (§io) 



1 — 2cut 



X fcj, fc Q ;TCr 
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)f (4-43) 
Mq (ki(x),k a (x);Tai^ 



The z-integrations and the sum over r can again be performed, resulting in 
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and 
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(4.45) 



4.3 The quark-to-photon fragmentation function 

In the previous sections we described how we can deal with identified particles in the final state that lead 
to non-IR-safe observables using phase-space slicing and the subtraction method. To restore IR safety 
one factorises the resulting left-over singularities into an experimentally determined fragmentation 
function, in our case the quark-to-photon fragmentation function. 

In Ref. [83] a method has been proposed in how to measure the quark-to-photon fragmentation 
function, i.e. the probability of a quark splitting into a quark and a photon, using the e + e~ — > 
njet + photon cross section. The method has been extended to next-to-leading order in QCD in 
Ref. [87]. The key feature of the proposed method is the democratic clustering of both hadrons and 
photons into jets, where one keeps track of the fraction of photonic energy in the jet. This treatment 
of the photon in the jet enhances the non-perturbative part of the quark-to-photon fragmentation 
function [25,83], which in turn can be measured in e + e _ annihilation. 

In Ref. [83] the quark-to-photon fragmentation function was theoretically defined using dimensional 
regularisation and one-cut-off slicing. To be able to use the results obtained in this way in our calcu- 
lation, we need to translate this definition to mass regularisation. We first summarize the results in 
dimensional regularisation. 

Since fragmentation is a long-distance process, it cannot be calculated entirely in perturbation 
theory. The fragmentation function D 9 _ >7 (^ 7 ) describes the probability of a quark fragmenting into 
a jet containing a photon carrying z 7 of the jet energy. Photons inside hadronic jets can have two 
origins: (a) The perturbatively calculable radiation of a photon off a quark, which contains a collinear 
divergence, described by a perturbative contribution dcr^ oll , dependent on the method used to regulate 
the collinear singularity; (b) the non-perturbative production of a photon in the hadronisation process 
of the quark into a hadronic jet, which is described by a bare non-perturbative fragmentation function 
Dq!^ 7 (z 7 ). Both processes refer to the infrared dynamics inside the quark jet, and can a priori not be 
separated from each other. 

Within dimensional regularisation and one-cutoff-slicing, the photon fragmentation contribution 
was studied in detail in Refs. [83,84]. Exploiting the universal factorisation of matrix elements and 
phase space in the collinear limit, one obtains for the cross section for the emission of one collinear 
photon with energy fraction z 7 above z cut off a quark q and invariant mass of the photon-quark pair 
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(4.46) 



where doo is a reduced cross section with the quark-photon system replaced by its parent quark. The 
same expression was obtained in (4.15) above for the hard-photon cut contribution, where z — 1 — z 1 . 

The infrared singularity present in this perturbative contribution must be balanced by a divergent 
piece in the bare fragmentation function, which contributes to the photon-emission cross section as 



dcr fr . 



;(Zcut) = d(T / 



dz 7 ^ 7 (z 7 ) 



(4.47) 



To make this cancellation explicit, one introduces a factorisation scale /x F into the bare fragmentation 
function, which then reads in dimensional regularisation 



n barc,DR, x n , „ x , a Q\ 1 / 4 ^ 
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(4.48) 



The factorised non-perturbative fragmentation function D q ^ J (z 1 , /x F ) is then finite, and can be de- 
termined from experimental data. Its variation with the scale /i F is described by the Altarelli-Parisi 
evolution equation, which reads to leading order 



dD q ^. 1 (z 1 ,fi F ) _ 



d ln /up 2-7T 
The fixed-order exact solution at O(a) is given by 



OtQ f u 2 \ 
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where I? 9 _j. 7 (z 7 , /iq) is the quark-to-photon fragmentation function at some initial scale (Xq. This 
function and the initial scale (io cannot be calculated and have to be determined from experimental 
data. The first determination of D q ^(z 7 fi ) was performed by the ALEPH collaboration [88] using 
the ansatz 
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with fitting parameters C and /xo- The fit to the photon-plus-one-jet rate [88] yielded 

it = 0.22GeV; C = -12.1 . 



(4.51) 



(4.52) 



Through the factorisation formula (4.48), a relation between bare fragmentation function and the 
collinear photon contribution is established. As a result, the sum of both contributions is finite, but 
depends on the slicing parameter s m i n : 



d^coii^cut) + dcr frag (z cut ) = da / dz 7 
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(4.53) 
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Inserting the solution (4.50) into (4.53), d<7^ oll + dcrf ra g becomes independent of the factorisation scale 
Mf- 

We can now determine the contribution of the quark-to-photon fragmentation function to 3-jet 
production. Since photons above z cu t are vetoed we have to subtract the photon-fragmentation con- 
tribution from the real corrections, with the hard-photon cut procedure included, and the virtual 
corrections, resulting in the IR-safe cross section 



J dcr IR - sato = J der virt + J dcr rcal (z cut ) - J dcr frag (z cut ). 



(4.54) 



Taking into account all quarks (and antiquarks) in the final state, the photon-fragmentation contribu- 
tion reads: 



f 1 bare 
d<7f ra g(;Z cu t) = y^d(TBorn / dz 7 Dg*^ 



- qi ^)- (4-55) 

Physically we can motivate this approach as follows. The hard photon cut removes collinear quark- 
photon pairs above a predefined photon energy fraction. This results in an incomplete cancellation of 
collinear singularities between the real photon radiation and the virtual QED-type corrections (where 
the hard photon cut is not acting, since the virtual photon is not observed). By subtracting dcrf rag 
from the O (a) corrections, we correct for the effect of the hard-photon cut and compensate for excess 
terms related to collinear photon emission. In this way we can define an IR-safe quantity even in the 
presence of the hard-photon cuts. 

In order to use the fragmentation function in our calculation with mass regularization we have to 
translate (4.48) to mass regularization. For consistency with previous sections, we use from now on 
the quark energy z — 1 — z 7 as collinear variable. 

Using results of Ref. [79], the collinear photon contribution in mass regularisation and one-cutoff 
slicing is obtained as 
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(4.56) 



Using this result and the independence of (4.53) on the regularisation scheme we find the bare frag- 
mentation function in mass regularization 
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where the finite terms ensure that the MS scheme factorised fragmentation function D q ^ 1 {z 11 /xf) is 
identical in the dimensionally regularised and the mass-regularised expressions. After inserting the 
ALEPH ansatz (4.51) for Z? 9 ^ 7 (z 7 , fip), we obtain 
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Integrating (4.58) over z results in 
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In the case of the two-cut-off phase-space-slicing method, subtracting d(7f r ag(z C ut); i-e. (4.55) with 
(4.59) from dcr^f (z cu t), (4.13), leads to 
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with dcr^"^ 1 from (4.5). This result consists of the original collinear contribution that cancels against 
the virtual corrections, and an additional, finite contribution, depending on the cut on the quark energy 
Zcut- All collinear singularities have cancelled. 

In the case of the subtraction method, we can use charge conservation 



Qi = - £(-i) (i+J) QiQj = -£(-i) (<+o) QiQa - ^(-i) (i+j) g i Q j 



(4.61) 
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in (4.55) to split d<jf lag (z cut ) into a contribution from final-state emitter and final-state spectator and 
one from final-state emitter and initial-state spectator 
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(4.62) 

In the case of a final-state emitter and final-state spectator, we can subtract the second part of (4.62) 
from the contribution of (4.35) to find 
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Analogously, in the case of a final-state emitter and initial-state spectator using (4.44), we are left with 
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Both (4.63) and (4.64) are finite and only depend on the value of z cut , but not on the mass regulators. 
4.4 Higher-order initial-state radiation 

In order to achieve an accuracy at the per-mille level, we also include effects stemming from higher-order 
ISR using the structure-function approach as described in Refs. [81,89]. The factorisation theorem 
states that the leading-logarithmic (LL) initial-state QED correction can be written as a convolution 
of the lowest-order cross section with structure functions according to 

J da LL = J dx,^ dx 2 T^(x u Q 2 )T^(x 2 ,Q 2 ) J da Boin (x 1 k 1 ,x 2 k 2 ), (4.65) 

where x\ and x 2 denote the fractions of the incoming momenta just before the hard scattering, Q 2 
is the typical scale at which the scattering occurs, and the structure functions up to 0(a 3 ) are given 
by [89] 
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T is the Gamma function, and 7e the Euler-Mascheroni constant. In the calculation at hand we use 
Q 2 = S . 

When we add (4.65) to the one-loop result, we have to subtract the lowest-order and one-loop 
contributions d<r LL1 already contained in (4.66) to avoid double counting. They read 
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where the one-loop structure functions are given by 
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5 Implementation 



The real and virtual corrections described in the two previous sections are implemented in a parton- 
level event generator. This program generates final states with two and three particles for the hadronic 
cross sections, and with three and four particles for event-shape distributions. It allows for arbitrary, 
infrared-safe cuts on the final-state particles to be applied. 

5.1 Event selection 

The infrared-safety requirement has one important implication that was up to now not necessarily 
realized in the experimental studies. In ISR events, where a photon is radiated close to the beam-pipe 
and not observed, the event-shape variables must be computed in the CM system of the observed 
hadrons. Only this transformation to the hadronic CM system ensures that two-jet-like configurations 
are correctly identified with the kinematic limit y — > of the event-shape distributions, as can be seen 
on the example of thrust: a partonic final state with a quark-antiquark pair and an unobserved photon 
yields two jets which are not back-to-back in the c+c~ CM frame. The reconstructed thrust axis in 
this frame is in the direction of the difference vector of the two jet momenta, resulting in T < 1 , even 
for an ideal two-particle final state. 

Unfortunately, most experimental studies of event-shape distributions at LEP computed the shape 
variables in the e + e~ CM frame. Modelling the ISR photonic corrections using standard parton-shower 
programs, these distributions were then corrected (by often large correction factors, shifting the two-jet 
peak from within the distribution back to the kinematical edge) bin-by-bin. It is therefore very difficult 
in practise to compare our results with data from LEP. Despite its importance this task goes beyond 
the present paper. 

In our implementation, we proceed as follows: 

1. For all final-state particles, an acceptance cut is applied on the polar angle 6i of the particle i 
with respect to the beam direction: |cos0j| < cos6> cut . Particles that do not pass this cut are 
unobserved and thus discarded, i.e. their momenta are set to zero. 

2. For the remaining final-state particles, the observed final-state invariant mass squared s' is com- 
puted. The event is rejected if s' < s cut . 

3. If the event is accepted, it is boosted to the CM frame of the observed final-state particles. 

4. To identify isolated-photon events, all observed final-state particles (including the photon) are 
clustered into jets using the Durham algorithm with E recombination and y cut = 0.002. After 
the clustering, the photon is inside one of the jets (or makes up one jet), carrying a fraction z of 
the jet energy. If z > z cut , the photon is considered isolated and the event is rejected. 

5. All remaining events are accepted. 

Once an event passes the above set of cuts, we proceed with the calculation of the event-shape 
variables, defined in Section 2, in the CM frame of the observed final-state particles. We impose an 
additional cut individually for each histogram such that the singularity in the two-jet region is avoided, 
typically as a lower cut-off for the respective observables (see Section 5.2). The cut affects only the 
first bin of the histogram and does not cause a distortion of the shape of the distribution. 

5.2 Input parameters 

For the numerical calculation, we use the following set of input parameters [90]: 



G M = 1.16637 x 10~ 5 GeV 
a s (M z ) = 0.1176, 
Mvf p = 80.403 GeV, 
M^ EP = 91.1876 GeV, 



—2 



a(0)= 1/137.03599911 



a Gfi = 1/132.49395637 



m c = 0.51099892 MeV, 



r^ EP = 2.141 GeV, 
r^ EP = 2.4952 GeV, 
m t = 171.0 GeV, 



Mh = 120 GeV. 



(5.1) 
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We employ the complex-mass scheme [65] , where a fixed width enters the resonant W- and Z-boson 
propagators, in contrast to the approach used at LEP to fit the W and Z resonances, where running 
widths are taken. Therefore, we have to convert the "on-shell" values of My EP and T\P P (V = W, Z), 
resulting from LEP, to the "pole values" denoted by My and Ty . The relation between the two sets is 
given by [91] 

My = ML EP / v /l + (rLEP/ M LEP)2 j Ty = T^/y/l + (T™ P /M™ P ) 2 , (5.2) 

leading to 

M w = 80.375 . . . GeV, T w = 2.140 . . . GeV, 

M z = 91.1535... GeV, T z = 2.4943 ... GeV. (5.3) 

The scale dependence of a s is determined according to the two-loop running. The number of active 
flavours at Mz is np — 5, resulting in A5 = 0.221 GeV. The scale dependence is matched to two-loop 
order at the top threshold [92]. 

We neglect effects due to quark mixing and set the CKM matrix to unity. Throughout this work, 
we parametrise the couplings appearing in LO in the scheme, i.e. we use ag^, whereas we fix the 
electromagnetic coupling appearing in the relative corrections by a — a(0). If not stated otherwise, 
we use the parameters 

cos e cut - 0.965, s cut = 0.81s, z cut = 0.9, y cut = 0.002 (5.4) 

in accordance with the event-selection criteria used in the ALEPH analysis [53] and employ the Durham 
jet algorithm together with the E recombination scheme for the reconstruction of isolated photons (see 
Section 2). 

As mentioned in Section 5.1, we implement a cut such that the singularity in the two-jet region is 
avoided. This cut requires the variables T, p, Bt, %, and C to be greater than 0.005, whereas I3 
and y cut for cr 3 _j et are required to be greater than 0.00005. 

5.3 Checks of the implementation 

The reliability of our results is ensured as follows: 

• UV finiteness is checked by varying the reference scale /x of dimensional regularisation and finding 
that our results are independent of this variation. 

• IR finiteness is verified through varying the infinitesimal photon mass m 7 in mass regularization 
and observing that the sum of the virtual corrections and the soft-photonic corrections in both 
the slicing and subtraction approach is invariant. In dimensional regularisation the independence 
of \x is checked as for UV divergences. 

• Mass singularities related to collinear photon emission or exchange are shown to cancel between 
the virtual and the subtraction endpoint contributions by varying the small masses of the external 
fermions. 

• Two completely independent calculations have been performed within our collaboration. We find 
complete agreement of the results for <7h a d, jet rates, and event-shape distributions at the level 
of the Monte Carlo integration error which typically is at the per-mille level. 

Furthermore, we compare the results obtained with phase-space slicing and the subtraction method, 
which are completely independent techniques. In Fig. 9 we show the mutual agreement of both tech- 
niques for the NLO EW results for <7had, in Fig. 10 for the full O(a) results for the three-jet rate 
with y cu t = 0.0006 at y/s = Mz, and in Fig. 11 for the full 0(a) results for the thrust distribution at 
y/s = 206 GeV. Note that Fig. 9 and Fig. 10 refer to the corrections to c+c~ — > qq and e + e~ — > qqg, 
respectively, i.e. to two independent implementations. We find that within integration errors the slicing 
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y/s = M z 


do-i-dcTBorn T<K1 


y/s = 133 GeV 


dcTi-dcTBorn TO/J 




dCTBorn 


dCTBorn L'°J 


<°d m /nb 
<a°d ak /nb 

<a L d°/nb 

N LO+1 ,o.L L/nb 


38.2845(15) 
37.8541(2) 
25.729(3) 
27.341(3) 


-1.1 

-32.8 
-28.6 


0.068858(2) 
0.068348(2) 
0.06269(2) 
0.06208(2) 


-0.7 
-9.0 
-9.8 



^=172 GeV ^'d/ZT'" t % ] = 206 GeV ^^2'"" t % ] 



^hadVnb 


0.0276993(8) 




0.0170486(5) 




<a°d ak /nb 


0.0276780(8) 


-0.1 


0.0170626(5) 


0.1 


<a L d°/nb 


0.024770(7) 


-10.6 


0.015197(4) 


-10.9 


NLO +h .o.L L/nb 


0.024633(7) 


-11.1 


0.015127(4) 


-11.3 



Vi = 500GeV d Xtr" [ % ] ^= 1000 GeV dg ^°" [%] 



<°d m /nb 


0.00241881(7) 




0.00059139(2) 




<a°d ak /nb 


0.00238722(7) 


-1.3 


0.00056838(2) 


-3.9 


<a L d°/nb 


0.0020665(7) 


-14.6 


0.0004856(2) 


-17.9 


NLO +1 ,o.L L/nb 


0.0020585(7) 


-14.9 


0.0004836(2) 


-18.2 



Table 1: Total hadronic cross section <7had {y/s) f° r LEP1 and LEP2 energies, and for yfs = 500 GeV and y/s = ITeV. 



results become independent of the cut-offs for 8 S < 10~ 3 and S c < 10~ 4 and fully agree with the results 
obtained using the subtraction method. For the sake of clarity, we show only curves for values of the 
slicing parameters that lie on the plateau in Fig. 11. For values of the slicing parameters outside the 
plateau, the behaviour follows the same pattern as in Figs. 9 and 10. 

It turns out that the subtraction method is more efficient in terms of run-time compared to phase- 
space slicing. To obtain results of comparable statistical quality, the phase-space slicing method 
required about one order of magnitude more events than the subtraction method. We therefore used 
the subtraction method to obtain the results presented in the following sections, both for Chad and jet 
rates and event-shape distributions. 

6 Numerical results 

The numerical parton- level event generator described above can be used to compute the 0{a) elec- 
troweak corrections to the total hadronic cross section, to event-shape distributions, and to the three-jet 
rate. 

6.1 Results for the total hadronic cross section 

The total hadronic cross section as a function of the CM energy and the corresponding NLO electroweak 
corrections have been shown in Ref. [24]. Here we list some numbers that have been used to extract 
S a> i and 5cr,>2,LL, as defined in (2.17) and (2.21), which enter the calculation of normalised event- 
shape distributions and jet rates. We use the event selection as described in Section 5.1 with the cut 
parameters given in (5.4). In the following, 'weak' refers to the electroweak NLO corrections without 
purely photonic corrections. 

Table 1 shows the Born contribution to Chad m the first row, the weak 0(a) contribution in the 
second row, the full 0(a) contribution in the third row, and the full 0(a)+ h.o. LL contribution in the 
fourth row for LEP1 and LEP2 energies, as well as for ^/s = 500 GeV and y/s = 1000 GeV. We show 
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the absolute results in nanobarn in the second and fourth columns and the relative corrections in per 
cent in the third and fifths columns. The numbers in parentheses give the uncertainties from Monte 
Carlo integration in the last digits of the predictions. For most energies, the electroweak corrections 
are sizable and negative, ranging between —30% at the Z peak and about —9% at energies above. The 
numerically largest contribution is always due to ISR. Above 120 GeV the magnitude of the corrections 
is increased due to LL resummation of ISR, whereas it is decreased in the resonance region. The virtual 
one-loop weak corrections (from fermionic and massive bosonic loops) are moderate of a few per cent 
and always negative for Mz < y/s < 1 TeV. 

6.2 Results for the event-shape distributions and jet rates 

In the following we present the results of our calculation for the three-jet rate as well as for event-shape 
distributions as described in Section 2. We show our findings for y/s = Mz as used at LEP1 and the 
selected LEP2 energies 172 GeV and 206 GeV. To stress the relevance of our work for future linear 
colliders, we also show results for \fs = 500 GeV. 

The precise size and shape of the corrections depend on the observable y in question. However, 
they share the common feature that qqj final states contribute only in the two-jet region, typically for 
small values of y. 

In a first step, we show our results for the distributions normalised to er for \/s = Mz in Fig. 12. 
The Born contribution is given by the A term of (2.18), while the full 0(a) corrections contain the 
tree-level qqj contribution <5 7 and the NLO electroweak contribution 5a of (2.18). The T, p, Bt, 
Bw, and C distributions are weighted by the respective variable y, evaluated at each bin centre. The 
relative corrections in the lower boxes are obtained by dividing the respective contributions to the 
corrections by the Born distribution given by the A term. We observe large negative corrections due to 
ISR, and moderate weak corrections in all distributions. The corrections are mainly constant for large 
y (note that we plot 1 — T instead of T), where the isolated-photon veto rejects all contributions from 
qqj final states. Near the two-jet limit, the contribution from qqj final states dominates the relative 
corrections. Moreover, it turns out that the electromagnetic corrections depend non-trivially on the 
event-selection cuts (see Section 6.3 for a more detailed discussion). We observe a significant decrease 
from the second bin to the first bin in all distributions, caused by the lower cut-off that we impose 
individually for all distributions. Since the cut-off acts both in the Born and the NLO contribution, we 
find a meaningful result for the relative corrections in the first bin. In the I3 distribution we clearly see 
the onset of the qqj final states for Y3 = 0.002. Since we always cluster photons with y < y cut = 0.002 
in the event selection (see Section 5.1), the contribution from qq-f final states is removed if I3 > 0.002 
and only plays a role for Y 3 < 0.002. 

In expanding the corrections according to (2.19), and retaining only terms up to LO in a s , we 
obtain the genuine electroweak corrections to normalised event-shape distributions, which we display at 
\/s = Mz in Fig. 13. The Born contribution is given by the A term of (2.19), while the 0{a) corrections 
are now given by Jew of (2.20). It can be seen very clearly that the large ISR corrections cancel 
between the event-shape distributions and the hadronic cross section when expanding the normalised 
distributions properly, resulting in electroweak corrections of a few per cent. Moreover, effects from 
ISR resummation are largely reduced as well, and the difference between 0(a) and 0(a) + h.o. LL is 
very small. The purely weak corrections are below the per-mille level. 

For the thrust distribution, the full O(a) corrections are almost constant around 0.5% for (1 — T) > 
0.05. The coefficient <5 7 starts to emerge for (1 — T) = 0.04 and contributes up to 2.6%. The full 0(a) 
corrections peak for (1 - T) = 0.02 at 2.5% and amount to 1.8% for (1 - T) = 0.01 in the first bin. 
The full 0(a) corrections are flat and around 0.5% for p > 0.05, B w > 0.05, C > 0.1, and Y 3 > 0.002. 
For the Bt distribution they are around 1% and almost flat for Bt > 0.05. The full O(a) corrections 
reach a maximum between 1% and 3% typically for small values of the event-shape variables and drop 
towards the first bin down to between 1.8% and —4.1%. Only for Bt we find also a maximum (of 4%) 
in the last bin. The LO qqj channel contributes only for small y and can amount up to 4%. 

In Fig. 14 we show our results for \/s = 172 GeV. The behaviour is similar as for \/s = Mz- but 
in the middle of the distributions a peaking structure emerges. It is located at 1 — T, p, £>w ~ 0.2, 
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Figure 12: The event-shape distributions normalised to <tq at y/s = Mz. 
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Figure 14: The event-shape distributions normalised to Uhad a t \/s = 172 GeV. 
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Figure 16: The event-shape distributions normalised to cr^d at y/s = 500 GeV. 
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Figure 17: The three-jet rate normalised to tr^ad at different CM energies. 



Bt ps 0.25, C pa 0.65, and I3 pa 0.1. We investigate this behaviour in detail in the next section. In all 
event-shape distributions the LO qqj contribution ranges between 3% and 8%. Outside the two-jet 
region and apart from the domain of the peaking structure, the full 0(a) corrections are flat and of 
the order of 5%. They peak near the onset of the qqj final states between 4% and 10%, and drop in 
the first bin down to between 1.5% and —10%. 

Results for y/s = 206 GeV, are displayed in Fig. 15. In all event-shape distributions the LO qqj 
contribution ranges between 4% and 9%. Outside the two-jet region and outside the domain where 
the peaking structure is located, the full O(a) corrections are flat between 0.1% and 2%, they peak 
near the onset of the qqj final states between 5% and 9%, and drop in the first bin down to between 
2% and —8%. The peaking structure is now situated at smaller values of y, it is less pronounced and, 
especially for B\y, Bt, C, and Y3, it extends over a larger range of y. Additionally, for large values of 
y, the weak contribution slightly increases. 

Finally, in Fig. 16 we show our results for y/s = 500 GeV. In the event-shape distributions the LO 
qqj contribution ranges between 2% and 8%. Outside the two-jet region, the full 0(a) corrections 
are flat between 2% and 3%, they peak near the onset of the qqj final states between 2% and 9%, 
and drop in the first bin down to between 2% and —6%. The weak corrections here contribute up to 
+1% in all observables for large values of y. The peaking structure as observed for y/s = 172 GeV and 
sfs = 206 GeV has completely disappeared. 

Figure 17 displays the corrections to the three-jet rate at various CM energies. As above, the 



37 



corrections are appropriately normalised to Chad- At y/s = Mz, the full 0(a) corrections to the three- 
jet rate are about 0.5% for y cut > 0.002. Because we use y cut — 0.002 in the event selection, ggj final 
states contribute only if y cut < 0.002. The LO qqj contribution amounts to 1%. For y cut < 0.002, the 
full O(o) corrections become negative, reaching —1.5% in the first bin. For y cut < 0.0005 the three-jet 
rate becomes larger than <7 n ad- This behaviour indicates the breakdown of the perturbative expansion 
in a s due to large logarithmic corrections proportional to log(j/ cut ) at all orders. Inclusion of higher- 
order QCD corrections, which are large and negative [17, 18] in this region, yields a ratio of cr3_j et to 
Chad less than unity for an extended range of log(j/ C ut)- At the higher CM energies, the corrections to 
the three-jet rate are larger than at y/s = Mz- For y cut > 0.002, they are almost constant and amount 
to about 4% at 172 GcV, and about 2% at 206 GeV and 500 GeV. In the region y cut < 0.002, we find 
a negative contribution of up to —5% for very small values of y cu t- 

6.3 Impact of the event-selection cuts on the event-shape distributions 

In the above results, we could clearly observe that the electroweak corrections to event-shape distribu- 
tions are not flat, but display peak structures. These structures are most pronounced at y/s = 172 GeV. 
They are discussed here for thrust as an example. In Fig. 14, we see that the relative corrections show a 
peaking structure for 1 — T s=s 0.2. To understand the origin of these structures, we extensively studied 
how variations of the event-selection cuts, especially the hard-photon cut, influence the event-shape 
distributions. 

We employ three different cuts in our calculation which depend on four parameters: 

1) A cut on the production angle 6i of all particles, such that only particles i with cos#j > cos# cut 
are used in the reconstruction of the event-shape variables. By default, we use the value cos 9 cut = 
0.965. 

2) A cut on the visible energy squared s' of the final state, such that only events with s' > s cut are 
accepted. By default, we use the value s cut = 0.81s. 

3) The maximum photon energy in a jet z, for which we require z < z cut , where the particles are 
clustered according to the Durham jet algorithm with parameter y cut . By default, we use the 
values z cut = 0.9 and y cut = 0.002. 

In Fig. 18 we show the full 0(a) corrections to the thrust distribution normalised to the Born 
contribution for y/s = 172 GeV, where the peak structures are most striking. We plot the results for 
three different values of a single cut parameter while we set the other three cut parameters to their 
default value. Going from left to right and top to bottom, we vary s cut , cos# cu t, z cu t, and y cu t- By 
changing s cu t we observe a change in normalisation of about 25%, while the shape of the distribution 
stays the same. Varying cos 6> cu t from larger to smaller values leads to a more and more pronounced 
peak at (1 — T) ~ 0.2. The corrections grow below the peak but are only slightly changed above. 
The different slopes for 1 — T < 0.03 result from the changing acceptance of ISR photons. Modifying 
z cu t has a dramatic effect on the peak. For z cut = 0.99 we find a very pronounced resonance for 
(1 — T) ~ 0.28. By reducing the cut, the resonance gets strongly suppressed and moves to smaller 
values of (1 — T). By increasing y C ut, we observe an enhancement of the peak, as well as a slight shift 
towards larger values of 1 — T. 

In Fig. 19 we study the change of the 0{a) corrections to the thrust distribution normalised to the 
Born contribution with z cut for y/s = Mz, 133 GcV, 206 GeV, and 500 GeV. Varying z cut for y/s = Mz 
leaves the distribution basically unchanged. By increasing z cut from 0.5 to 0.99 for y/s — 133 GeV we 
find a growth of the peak to almost 100%. We also observe that the peak moves from 1 — T ss 0.28 
for z cut = 0.9 to 1 - T w 0.25 for z cut = 0.99. For y/s — 206 GeV we see the same features as for 
y/s = 172 GeV, but now with the peak at (1 - T) ~ 0.19 for large z cut . Finally, for y/s = 500 GcV 
a peaking structure seems to appear at (1 — T) ~ 0.03, while varying z cut basically changes only the 
normalisation. 

Through analysing events at the level of the Monte Carlo generator, we find that the enhancement 
in the region of the peaking structure always stems from gggj final states, where a soft gluon is 
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Figure 18: Dependence of the thrust distribution on different values of the phase-space cuts at y/s = 172 GeV. 



clustered with a hard photon. Increasing cos 9 cut leads to a logarithmic enhancement of collinear ISR 
photons, increasing z cut generally results in a larger acceptance of photons inside jets, and increasing 
y cut causes more photons to be clustered together with other partons, resulting in less events with 
isolated photons being removed. We can therefore conclude that the peaking structure results from 
the ISR photon contribution, where a soft gluon is clustered together with the photon. 

More precisely, the peaking structure can be explained by the radiative-return phenomenon. Since 
we do not remove all energetic photons, it is possible that a hard photon and a soft gluon are clustered 
together, such that the energy fraction of the photon in the jet does not exceed z cu t and the invariant 
mass of the quark-antiquark-gluon system s q q g is equal to the mass of the Z boson. Such a configuration 
leads on the one hand to an enhancement due to radiative return but also to a logarithmic enhancement 
due to the soft gluon. 

In order to analyse this effect further, we consider events where the photon and the soft gluon are 
clustered together, such that we have a three-particle final state that consists of a quark, an antiquark, 
and a photonic jet. Assume that the quark, antiquark, and the photonic jet have the three-momenta 
PqiPqiPi an d the energies E q ,Eq,E 1 , respectively. We use energy-momentum conservation 

Pq+Pq+Pj = 0, 
E q + Eg + E y = yft. (6.1) 

and demand that the invariant mass of the quark-ant iquark pair is equal to Mz, such that 

Sqq = {E q + E q f - (p q + p q f = (E g + E q f - (p 7 ) 2 = {E q + E q f -E% = Ml (6.2) 
Energy conservation (6.1) and the mass-shell condition (6.2) imply 

s + M% s-MS , s 
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Figure 19: Dependence of the thrust distribution on different values of the cut z cu t at y/s = Mz, 133 GeV, 206 GeV, 
and 500 GeV. 



It can be shown that in a three-jet configuration with massless partons, thrust is always determined 
by the energy of the most energetic particle B max [27], i.e. 

T= 2E^ 

If we now assume that the quark and the antiquark carry the same energy, we can calculate the energies 
of all three jets in the final state at different energies: 

E q (133 GeV) = Eg (133 GeV) 49 GeV, E 7 (133 GeV) « 35 GeV, 

E q (172 GeV) = Eg (172 GeV) « 55 GeV, £ 7 (172GeV)« 62 GeV, 

E q (206 GeV) = E q (206 GeV) w 61 GeV, E y (206 GeV) w 84 GeV, 

S 9 (500 GeV) = E q (500 GeV) w 129 GeV, fJ 7 (500 GeV) w 242 GeV. (6.5) 

This leads to the following thrust values where the radiative-return phenomena should appear: 

(1 -lRR)(Vi= 133 GeV) w 0.27, 
(1 - T RR ){yG = 172 GeV) « 0.28, 
(1 - Trr)(v^ = 206 GeV) ps 0.19, 

(1 - 7r, r )(v7 = 500 GeV) « 0.03. (6.6) 

These values coincide perfectly with the peaks in Figs. 18 and 19. Relaxing the assumption that the 
quark and the antiquark carry the same energy only results in a broadening of the peaking structure. 
Varying the value of z cu t leads to different energies of the photonic jet and therefore changes the 
allowed energies in the above analysis. For decreasing values of z C ut, we therefore either only observe 
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Figure 20: Electroweak corrections to the thrust distribution at different CM energies. The interference contribution 
between electroweak and QCD diagrams for the four-quark final state is scaled by a factor 1000. 



the tail of the peak or cut it away completely which effectively looks like a shift of the position or the 
disappearance of the peak. For y/s — 133 GeV we observe the tail of the peak for E q < E q ^ ca k such 
that for decreasing z cu t the peak seems to move to larger values of (1 — T). For s/s = 172 GeV and 
y/s = 206 GeV we observe the tail of the peak for E~ t > -E 7 . pGa k such that for decreasing z cut the peak 
seems to move to smaller values of (1 — T). 

The study in this section clearly illustrates the non-trivial effect of realistic photon isolation criteria 
to the electroweak corrections to jet observables. The accidental clustering of a soft gluon with a hard 
photon results in a photon jet with a photon energy fraction below the rejection cut. In these events, 
the distribution of the final-state jets, and their reconstructed pair invariant masses do no longer 
reflect the underlying parton dynamics. A similar misreconstruction could also happen for electroweak 
corrections to final states involving jets at hadron colliders, and clearly deserves further study. 

6.4 Contribution from four-quark final states 

At 0(a 3 a s ), event-shape distributions and jet cross sections receive a contribution from the process 
e + e~ — > qqqq through the interference of QCD and electroweak amplitudes [see (3.23)]. Compared to 
other contributions at this order, this four-quark interference contribution is very small. Its typical 
magnitude amounts to about one per mille of the electroweak correction, and is thus within the inte- 
gration error of the four-particle contribution. To illustrate the smallness of this contribution, Fig. 20 
compares the four-quark contribution (scaled by a factor 1000) to the total electroweak correction to 
the normalised thrust distribution at different CM energies. The relative magnitude of the four-quark 
contribution is always at the per-mille level. 



41 



7 Conclusions and outlook 



In this paper, we have derived the NLO electroweak corrections to three-jet production and event- 
shape distributions in e+c~ annihilation. At this order, contributions arise from virtual corrections 
from weak gauge bosons (which were evaluated in the complex-mass scheme to take proper account 
of the gauge-boson widths), from fermionic loops, from real and virtual photonic corrections and from 
interferences between electroweak and QCD diagrams for four-quark final states. 

Our calculation is one of the first to address electroweak corrections to jet production observables. 
For this type of observables, one has to take proper account of the experimental event-selection cuts, 
which are aiming to reject events containing isolated photons. An infrared-safe definition of isolated 
photons must permit some amount of hadronic energy around the photon; for jet observables, this can 
be realized by a democratic clustering procedure, used by the LEP experiments. In this approach, 
the photon is clustered like any other hadron by the jet algorithm, resulting in a photon jet. If the 
photon carries more than a predefined large fraction of the jet energy, it is called isolated and the event 
is discarded. In our calculation, we have implemented this isolated-photon-veto procedure. Since it 
involves cuts on a specific, identified particle in the final state, the resulting observable is no longer 
collinear-safe, and the left-over collinear singularity is compensated by a further contribution from the 
quark-to-photon fragmentation function. We have documented this part of the calculation in detail in 
different schemes. 

The NLO electroweak corrections to absolute cross sections and event-shape distributions turn out 
to be numerically substantial: for example for thrust at y/s — Mz, they amount to a correction of 
—32%, which is largely dominated by initial-state radiation. Beyond the NLO electroweak corrections, 
we also included higher-order leading logarithmic corrections, which are sizable. Their inclusion at 
yfs — Mz results in a total correction of —28%. Normalizing these results to the total hadronic 
cross section (as is done in the experimental measurement), corrected to the same order, only very 
moderate corrections remain in the normalized jet cross sections and event distributions, and practically 
no difference is observed between the fixed-order NLO electroweak results and the results including 
higher-order logarithmic corrections. 

At LEP1, we find that NLO electroweak effects to event-shape distributions amount to a correction 
of about one to two per cent. The corrections are not uniform over the kinematical range, but tend 
to increase towards the two-jet region, where the isolated-photon veto becomes less efficient. The 
corrections to the three-jet rate are below one per cent. Purely weak contributions form a gauge- 
invariant subset of the electroweak corrections. At LEP1, these corrections are below the per-mille 
level. 

At LEP2 energies, the NLO electroweak corrections to event-shape distributions and to the three-jet 
rate are typically at the level of two to eight per cent. The largest contribution comes again from the 
photonic corrections. These are very sensitive to the precise photon isolation cuts applied to select the 
events, and are not uniform over the range of event-shape variables. The NLO electroweak event-shape 
distributions display peaks at LEP2 energies. These peaks are due to a remant of the radiative-return 
phenomenon, which is not fully suppressed by the photon isolation cuts. The position and energy 
dependence of these peaks can be explained quantitatively. 

Event-shape and jet cross section data from LEP1 and LEP2 have been corrected for photonic 
radiation effects using multi-purpose leading-logarithmic event generator programs. To compare our 
results with the experimental data would first require to undo these corrections. A further compli- 
cation in the comparison with data arises from the fact that event-shape distributions at LEP2 were 
determined in the e + e~ centre-of-mass frame for all events, including initial-state-radiation events. 
Most of the event-shape variables are not boost invariant, and should thus be reconstructed in the 
centre-of-mass frame of the observed hadrons. If reconstructed in the e + e~ centre-of-mass frame, ideal 
two-jet events with initial-state radiation will not be placed at the kinematical edge of the distribution, 
thus violating infrared-safety criteria. Within a perturbative calculation, it is not possible to apply 
the event reconstruction in the e + e _ centre-of-mass frame. The purely weak corrections at LEP1 and 
LEP2 were previously not accounted for in the interpretation of event-shape and jet cross section data. 
Our study shows that they are at the level of one per mille or below for appropriately normalized 
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distributions. At the current level of experimental and theoretical precision, they are thus not yet 
relevant to precision QCD studies of LEP data. 

While the magnitude of electroweak corrections decreases towards higher LEP2 energies, we observe 
them to increase again when going to even higher energies, corresponding to a future linear collider. In 
part, this increase comes from the fact that purely weak corrections (which were negligible throughout 
the LEP2 energy range) become sizable at high energies. At y/s = 500 GeV, NLO electroweak cor- 
rections to event-shape distributions and jet cross sections amount to two to four per cent, and have 
thus a potentially sizable impact on precision QCD studies at a future linear collider. The purely weak 
corrections reach up to one per cent at this energy. Most importantly, our findings on the interplay 
of photon isolation and event-selection cuts, and on the appropriate frame for the reconstruction of 
initial-state radiation events will help to optimize precision QCD studies at future high-energy c + c~ 
colliders from event shapes and jet cross sections. 
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